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Abstract 

In various statistical-mechanical models the introduction of a metric onto the 
space of parameters (e.g. the temperature variable, (3, and the external field variable, 
h, in the case of spin models) gives an alternative perspective on the phase structure. 
For the one-dimensional Ising model the scalar curvature, 7Z, of this metric can 
be calculated explicitly in the thermodynamic limit and is found to be 1Z = 1 + 



cosh(/i)/y sinh 2 (/i) + exp(— 4(3). This is positive definite and, for physical fields and 
temperatures, diverges only at the zero-temperature, zero-field "critical point" of 
the model. In this note we calculate 7Z for the one-dimensional (/-state Potts model 
finding an expression of the form 1Z = A(q,(3,h) + B{q, (3, h)/^ r](q, (3, h) , where 
rj(q,(3, h) is the Potts analogue of sinh 2 (/i) + exp(— 4(3). This is no longer positive 
definite, but once again it diverges only at the critical point in the space of real 
parameters. We remark, however, that a naive analytic continuation to complex 
field reveals a further divergence in the Ising and Potts curvatures at the Lee- Yang 
edge. 



1 Introduction 



The notion of a distance between configurations in statistical-mechanical and lattice-field 
models has been discussed by several authors 0, §, f|, || and the geometry that 
this endows on the manifold, Ai, of the parameters characterising the models explored. 
This distance is measured using the Fisher-Rao metric J/J, which is calculated from the 
Fisher information matrix for the system of interest. When we consider a spin model in 
field, M. is a two-dimensional manifold parametrised by (6 l , 6 2 ) = (/3, h). In this case, the 
components of the Fisher-Rao metric take the particularly simple form 

Gij = -didjf, (1) 

where / is the reduced free energy per site and d{ = d/dO' 1 . 

It has been suggested that in such a geometrization of statistical mechanics the scalar 
curvature, TZ, of M. plays a central role |2|, [|. For a spin model in field, with the metric 
given in equ.(|l|), TZ may be calculated succinctly as 



TZ 



2G 2 



djf d p d h f dlf 
d$f dp h f dpdlf 
dp h f dpdlf dlf 



(2) 



where G = det(Gy). Since the only scale present near criticality for a model displaying a 
higher-order transition is the correlation length, £, it has been hypothesised on dimensional 
grounds that TZ ~ £ d , where d is the dimensionality of the system |], 0, If we assume 
that hyperscaling holds, vd — 2 — a, this leads to 

n ~ |e|( 2 - Q )/-. (3) 

To test the behaviour of TZ, one requires models which are solvable in field in order to 
obtain analytic expressions, and these are rather thin on the ground. Indeed, TZ has been 
calculated for the mean-field and Bethe-lattice Ising models and the above scaling 
behaviour verified. It has also been calculated for the one-dimensional Ising model 
where it takes the remarkably simple form, 

^ , cosh h . . . 

ft Ising = 1 + (4) 



In this case TZ is positive definite and diverges only at the zero-temperature, zero-field 
"critical point" of the model. The correlation length is given by 

T 1 = -ln(tanh(/3)), (5) 

so that £ ~ exp(2/3) near criticality, and (||) holds there with a = 1, v = 1 as expectecf]. A 
second noteworthy feature of the Bethe-lattice and ID Ising models concerns the metric 
associated with ([[]). This metric is diagonalized using the corresponding renormalization 
group invariant. 

lr The Bethe lattice model also satisfies the postulated scaling, although there are some subtleties 
coming from the exponent a being zero fl . 



Given the shortage of explicitly calculable examples, any further additions to the list 
would be very worthwhile in order to see which features in the models are generic and 
which are particular to the models concerned. In this paper we discuss another example, 
the ID g-state Potts model, where an expression for 1Z may be obtained in a very similar 
manner to the Ising model (which is the q = 2 case of the Potts model). We compare 
the properties of the curvature, TZ, as well as the metric in the Potts model and the Ising 
model, highlighting both their structural similarities and differences in detail. 



2 The ID Potts Model 



The partition function for the ID g-state Potts model is given by 



Z N (y,z) =J2 ex P 
M 



(6) 



where the spins, Oj G {1, . . . , q}, are defined on the sites, j G {1, . . . N}, of the lattice and 
where we have defined y = exp(/3) and z = exp(h) for later calculational convenience. 
The model may be solved by transfer matrix methods ||, just as the ID Ising model. 
For general q the full transfer matrix T(y, z) may be written as q — 2 diagonal elements, 
(y — l){yz)~ l / q , and a 2 x 2 factor t(y, z): 

z^y/q^T 



yz 



(yzy/tKz 1 / 2 ^^! y + q-2 



(7) 



The partition function is Z N (y,z) = TrT(y, z) N and the eigenvalues of T(y,z) are 
A ,Ai,...,A,_i, where 



A 
Ax 



y(l + z) + q - 2 ± J(y(l - z) + q - 2f + (q - l)Az 



(8) 



and A2, • • • A g _i = (y — l)(yz)~ l l q . The reduced free energy per site in the thermodynamic 
limit, iV — > 00, is thus given by / = — lnAo- 

It is straightforward to use this expression for the free energy in equ.(|2|) to obtain the 
curvature, 1Z. In the current notation we re-derive the expression for the Ising model | as 

9(1 + 2) (9) 



n 



Ising 



1 + 



\ly 2 — 2y' 2 z + y' 2 z' 2 + Az 
The expression for general q is similar in form to this, and is 

B(q,y,z) 



ftpotts = A(q, y, z) + 



(10) 



where the coefficients may be further broken down as A(q, y, z) = a(q, y, z)/ r y(q, y, z) 2 and 
B(q,y,z) = P(q,y, z)/j(q,y, z) 2 and are smooth functions of y and z and do not diverge 
for finite (physical) temperature or field. Furthermore 

V(q,y,z) = [y{l -z)+q-2] 2 + (q- 1)4*. (11) 

2 There is a factor of two difference in the definitions of /3, h between the Ising and Potts notations 
coming from the different spin definitions. 



The expressions for a(q,y,z), (3(q,y,z) and j(q,y,z) are very lengthy for general q (al- 
though still easily obtained) so, for the sake of compactness, we write down only those 
for q = 3, which already display the generic features. We have 

a(3, y,z) = --[-6-12z+ 116 z 2 + 80 z 3 - 16 z 4 

+y(-44 -68z + 400 z 2 + 552 z 3 - 192 z 4 ) 
+y 2 (-113 - 64 z + 138 z 2 + 476 z 3 - 1004 z 4 ) 
+y 3 (-178 + 30 z + 384 z 2 - 1884 z 3 - 2896 z 4 + 8 z 5 ) 
+y 4 (-188 + 458 z - 222 z 2 - 276 z 3 - 3400 z 4 + 64 z 5 ) 
+y 5 (-106 + 592 ^ - 474 z 2 - 68 z 3 + 1232 z 4 + 120 z 5 ) 
+y 6 (-17 + 294 z - 368 z 2 + 250 z 3 - 127 z 4 - 32 z 5 ) 

+y 7 (4 + 66 z - 136 z 2 + 60 z 3 + 2 z 5 + 4 z 4 )] , (12) 



and 



/5(3, y, z) = ~ [-30 + 12 z + 468 z 2 + 768 z 3 + 240 z A 

+y(-82 - 82 z + 664 z 2 + 3604 z 3 + 1744 z 4 - 16 z 5 ) 

+y 2 (-171 + 12 z + 1366 z 2 + 2292 z 3 + 1764 z 4 - 160 z 5 ) 

+y 3 (-261 + 213 z + 1710 z 2 + 846 z 3 - 7584 z 4 - 756 z 5 ) 

+y 4 (-278 + 386 z + 400 z 2 - 3804 z 3 - 8804 z 4 - 2488 z 5 + 8 z 6 ) 

+y 5 (-254 + 570 z - 1094 z 2 - 5654 z 3 - 1584 z 4 - 3712 z 5 + 64 z 6 ) 

+y 6 (-165 + 504 z - 1494 z 2 - 1008 z 3 + 891 z 4 + 1152 z 5 + 120 z 6 ) 

+y 7 (-51 + 263 z - 866 z 2 + 986 z 3 - 203 z 4 - 97 z 5 - 32 z 6 ) 

+y 8 (-4 + 66 z - 182 z 2 + 188 z 3 - 72 z 4 + 2 z 5 + 2 z 6 )] , (13) 

for the numerators, and 

7(3, y,z) = — 3 — 4z — 2z 2 — 4y — 16yz — 16yz 2 

- 7y 2 + 2y 2 z-31y 2 z 2 -4y 3 + 4y 3 z 2 , (14) 

and 

77(3, y,z) = l + 8z + 2y-2yz + y 2 -2y 2 z + y 2 z 2 , (15) 

for the terms in the denominators. 

In zero-field (z — 1) the expression for 72. is much more compact and is written for 
general q as 

(y + g - l)(4y 2 + (g - 2)y - (g - 2)(q - 1)) 

(g-l)(2y + g-2)2 " ^ 

We see that as y ranges from 1 to 00, 72 ranges from (4 — q)/ (q — 1) to 00. In particular, the 
sign of the y = 1 (@ — 0) limit of 72 changes at q — 4, although the general morphology 
of 72 as a function of y and z remains the same for all q > 2 as we see below. 



The correlation length for the one-dimensional Potts model is defined in a similar 
manner to that of the Ising model, 




r^-Mn (it) 



so £ ~ y for z — 1, y — > oo. We thus retrieve the expected scaling of TZ for the one- 
dimensional Potts model from equ.(|3]), namely TZ ~ y as y — > oo. The exponents, as for 
the Ising model, are a — 1, u — 1. 

The general features of 7?. at non-zero temperature and field are perhaps easiest seen 
in a contour plot as a function of ?/ and z. In Fig.l we show the Ising (q = 2) case which 
has certain non-generic features. The ±/i symmetry of the Ising model is manifest as a 
2 — > 1/z symmetry in the plot of TZ. In addition, one can see that TZ is positive for all y 
and z. The maximum of TZ for a given y value lies along the zero field line at z — 1. 

In Fig. 2, TZ is plotted for the 3-state Potts model, using the expressions given above 
for A(3,y,z) and B(3,y,z). We see that there is no longer a z — > 1/z symmetry and 
that 7£ is not positive definite. This behaviour is typical of the q > 3 case, as can also 
be seen from the q = 10 results plotted in Fig.3. It is also clear from Fig. 3. that the 
z — l,y — 1 limit is negative unlike the 3-state model, as indicated by equ.([Lj]). When 
q > 2 the maximum of TZ no longer lies on the zero-field line though its locus is still simply 
determined, as we discuss in the next section. 



3 Remarks on Co-ordinates and Geodesies 

We have concentrated so far on the curvature TZ, which is of course a geometric invariant 
and independent of the particular co-ordinate scheme we use in the calculation. The 
properties of the particular metrics and line elements used in the calculation are also of 
some interest ||. The metric for the ID Ising model calculated from equ.(|l|) gives the 
line element 



ds 2 



e -4/3 



(sinh 2 h + e- 4 0) 3/2 



4e 4/3 cosh h + 8 sinh /ifcosh h + \ sinh h + e 4/3 ) 9 

■ df3 2 



2 l, , —4/^2 



(cosh h + Y srnn h + e 



+ 4 sinh h d(3dh + cosh h dh 



Although this is clearly non-diagonal, the choice of (3 and h as co-ordinates has the 
advantage of giving the simplified expression, equ.(^), for TZ since they appear as linear 
multipliers in the Hamiltonian. 

The metric can be diagonalized by inspection in this case, and this preferred choice 
of co-ordinates turns out to be rather interesting from the physical point of view. If one 
considers a decimation type renormalization transformation on a ID Ising chain (or ring) 
where every other vertex is decimated and the length scale suitably adjusted, then one 
can exactly determine the relations between the renormalized parameters j3', h' and the 
originals (3, h from 

Z N/2 (f3',h>)=A N Z N (f3,h) (19) 

where is the decimated partition function and ZN(j3,h) is the original one. In 

addition to the parameter transformation there is also an unimportant overall scaling 



factor, A N . This transformation has an invariant p = e 2/3 sinh h, which is in effect encoding 
the invariance of the magnetization under the transformation since 



M 



VTT7' 

If the co-ordinates are transformed from /3, h to /3, p then, using 

dp = 2e 2/3 sinh hd(3 + e w cosh hdh, 
one finds the diagonalized expression 



ds 2 



1 



l + p 2 )(e^ + p 2 ) 



4 e^d(5 2 



+ 



dp 2 



(20) 



(21) 



(22) 



The use of the renormalization group invariant, or alternatively the magnetization, 
also diagonalizes the metric for the Ising model on a Bethe lattice. It is therefore natural 
to ask whether this feature persists in the Potts models considered here. Although the 
Potts metric and line elements are rather more complicated than those for the Ising model, 
the same general structure is once again apparent 



ds' 



X 



C(q, y, z) + D(q, y, z)Jij(q, y, z 



■d(3 2 + A[yz{z - l)}dpdh 



r](q,y,zy 3 /V ' [ \(q,y,z) 
+ z[y(l + z) + q - 2}dh 2 

where rj(q,y,z) is defined in equ. (|Tl"D and 

\(q, y,z)=y(l + z) + q-2 + J(y(l - z) + q - 2) 2 + (q - l)Az 



(23) 



(24) 



is proportional to Ao, the dominant eigenvalue. The functions C(q, y, z) and D(q, y, z) are 
both easily calculable, but rather long, and are not reproduced here. 
If we carry out a decimation on the ID Potts model 



then the invariant is given by 



J N/2 



P 



(y',z')=A»Z N ( y ,z) 



[y(l-z) + q-2] 



(25) 



(26) 



which is again related to the magnetization. Following the Ising procedure, we can change 
variables /3, h to /3, p using 



1 fz 2 x/i 



(27) 



and find that the metric is, indeed, diagonalized. 

In the Ising model it can be shown that the line p = (i.e. h = or z = 1) is a 
geodesic of the metric ([!]). In Fig.l this is the line running along the ridge of 1Z. For 



q 7^ 2 Potts models the local maximum in 1Z no longer lies at z = 1, as is clear from Fig. 2 
and Fig.3, however an analysis of the geodesic equations 



^+T%VPvt + 2Tl h VPV h + T p hh V h V h = \{s)V p 
dV h 

-J- + T^ p V p V p + 2T h ph V p V h + V h hh V h V h = X(s)V h (28) 

where s parameterizes the flow lines, = d/3/ds, V = dh/ds, the T are the various 
Christoffel symbols for the metric of equ. ([!]), and A(s) allows for non-affine parameters, 
shows that the line z = 1 is still a geodesic. To prove this consider a vector field with a 
flow line along z = 1 {h = 0). This flow line has V h = so equations (|28| ) become 



ds 



= 0. (30) 

h=0 K 1 



The first of these equations is a second order ordinary differential equation for j3(s) which 
always has a solution. The second requires (r^| h=0 = and it can easily be shown, 
using the metric (23) , that this is indeed the case. Hence the line z = 1 is a geodesic of 
the metric (^3[) for any value of q. 

A different change of variable in the Potts models brings one back to something very 
similar to the Ising picture. If we transform z to 

w = 31 

y+q-2 1 ; 

and leave y untouched, we find that the local maximum in 1Z lies on the line w — 1. 
We show 1Z plotted against w in Fig.4 for the 3-state Potts model for comparison with 
Fig.2. The picture is similar for higher q. This transformation does not greatly simplify 
cu(q, y, z), (3(q, y, z) or 7(5, y, z) and we do not reproduce the expressions here. 

In summary, although the choice of /3,h as parameters is natural in any calculation 
of 7Z, the metric is diagonalized for both the ID Ising and Potts models when the renor- 
malization group invariant, or equivalently the magnetization, is used as a co-ordinate 
instead of h. The line z = 1 is a geodesic of the metric ([I]) for any q. 



4 A Lee- Yang Divergence 

Some years ago Lee and Yang addressed the question of how the singularities associ- 
ated with field-driven phase transitions in Ising-like spin models on lattices arose in the 
thermodynamic limit. This was later extended by various authors to other models and 
to temperature-driven transitions jlO], |TTJ . Lee and Yang observed that the zeroes of the 
partition function for a spin model in a complex external field on a finite lattice would give 
rise to singularities in the free energy. In the thermodynamic limit these complex zeroes 
move in to pinch the real axis, signalling the the onset of a physical phase transition. 
Typically, the loci of zeroes are lines in the complex field or temperature plane and when 
the endpoints of such lines occur at non-physical (i.e. complex) external field values they 



can be considered as ordinary critical points with an associated edge critical exponent, 
usually dubbed the Lee- Yang edge exponent 

The Lee- Yang zeroes for the one-dimensional Potts model on a periodic chain with N 
sites are given by the solutions || [T^] of 



Z N = (X 1 ) N + {X o ) N = & Al = exp(^)A (32) 

where Aq,i are the eigenvalues given in equ. (|8|) and n is odd. In the thermodynamic limit 
the locus of zeroes is determined by |Ao| = |Ai| or 

r){q, y, z) = [y(l - z) + q - 2] 2 + (q - 1)4* = (33) 

which can be satisfied for complex (in the q = 2 Ising case, purely imaginary) values of h. 

Unlike the Ising model case, the Lee- Yang zeroes for the ID q ^ 2 Potts models do not 
lie on the unit circle in the complex z plane. However, precisely the change of co-ordinates 
we used in discussing the Potts metrics in the previous section, w = yz/(y + q — 2), places 
them on the unit circle in the complex w plane. If we associate a "field" with w via 
w = exp(h), then purely imaginary values of h give the zeroes, as in the Ising case?]. 

From these considerations, it is clear that TZ will also diverge as the locus of zeroes 
is approached for both Ising and Potts models if we allow ourselves the liberty of an 
analytic continuation of the field to complex h values once TZ has been calculated, since 
TZ = A + B/y/rj and A,B are finite as rj — > 0. The presence of the square root means 
that the divergence is characterised by an exponent a = — 1/2 which is the Lee- Yang edge 
exponent for the one-dimensional Potts (and Ising) model |T(| . 



The status of these observations are a little unclear to us, since the calculation of TZ 
has assumed a real metric geometry throughout and such an arbitrary continuation in the 
final expression might be rather dangerous. One is on slightly firmer ground with the Ising 
model, since in that case the required continuation is to purely imaginary fields which 
correspond simply to a change in the signature of the metric in /3, p co-ordinates. The use 
of the w co-ordinate in the Potts case does, however, suggest a similar interpretation since 
there the zeroes occur for imaginary values of h. With these caveats, it is nonetheless 
interesting that the Lee- Yang edge transition is still visible as a divergence in TZ. 



5 Conclusions 

We have seen that the scalar curvature, TZ, of the Fisher-Rao metric may be obtained for 
the ID g-state Potts model in a very similar fashion to the ID Ising model since the free 
energy can be calculated in field in both cases using transfer matrix methods. 

Although TZ for the g-state Potts model cannot be expressed as succinctly as for the 
Ising model it still has the same general form, TZ = A + B / ^/rj, and diverges only at 
the zero temperature critical point of the Potts model for physical temperature and field 
values, i.e. y — 1 . . . oo and z — . . . oo. We wrote down TZ explicitly here for q = 3, but 
noted it was a simple matter to calculate it for arbitrary q. We observed that there were 
some features of the Ising model TZ which did not persist for general q, as was clear from 

3 This is not, of course, the field h appearing in the Hamiltonian. 



contour plots. In particular, TZ was no longer positive definite and the z — * 1/z symmetry 
of the Ising model was no longer present. 

The choice of co-ordinate scheme used in the calculation of TZ was also discussed. We 
noted that although /3, h were a natural choice from the point of view of simplifying this 
calculation, the use of the renormalization group invariant in place of h diagonalized the 
metric for both the Ising and Potts models. A different choice of coordinates, involving a 
rescaling of z was also employed to identify the local maximum in TZ. For any value of q 
the line z — 1 (h — 0) is a geodesic of the Fisher-Rao metric. Finally, we observed that if 
complex field values are permitted there is a divergence in TZ at the Lee- Yang edge with 
an exponent a = —1/2. 

We have confined our discussion here to physical values of q, i.e. 2,3,.. ., in the Potts 
models. Since q appears as a parameter in the transfer matrix solution there is in principle 
no barrier to extending it to general non- integer q and also to q < 2. The q — > 1 limit of 
the Potts model is related to percolation, for example, so the behaviour of TZ might be 
of interest in this limit too. It would also be interesting to calculate TZ for models with 
genuine transitions at finite (3, an exercise that has so far been confined to mean-field or 
mean-field-like spin models. Possibilities which suggest themselves in this context are the 
Ising model on planar random graphs [T3] and the spherical model. 
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Figure 1: A plot of 1Z for the Ising model for y — 1 . . . 10, z — . . . 5. The positivity of 
TZ and the expected z — > \ jz symmetry are both apparent. 



Figure 2: A plot of 1Z for the 3-state Potts model for y — 1 ... 10, z — ... 5. 72. is 
no longer positive definite for physical values of y, z and there is no z — > \ jz symmetry. 
In addition the maximum of 1Z does not lie at z — 1, though this is perhaps clearer for 
higher q. 



10 



Figure 3: A plot of 1Z for the 10-state Potts model for y — 1 . . . 10, z — . . . 5. In this 
case it is quite clear that the maximum of 1Z does not lie at z — 1. 




Figure 4: A plot of 1Z for the 3-state Potts model for y = 1 ... 10, w = ... 3. Note that 
the maximum of 1Z lies at w — 1 



